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^Comparisons  have  been  made  between  a perfect  gas  laminar  viscous  shock 
layer  method  and  experimental  data  from  the  AEDC/vKF  Hypersonic  Tunnels.  The 
test  model  was  the  Rockwell  International  139  Shuttle  Orbiter,  and  the  data  i 

taken  included  wall  pressure  distributions,  heat-transfer  distributions, 

stagnation  temperature  and  pitot  pressure  profiles  through  the  shock  layer.  , 

Comparisons  were  made  at  Mach  numbers  of  8,  11  and  14,  with  Reynolds  numbers 

ranging  from  0-.&2-X  10-  to  1.0  per  foot  and  angles  of  attack  from  0 to 

45  degrees.  Comparisons  showed  good  to  excellent  agreement  except  at  very 

low  angles  of  attack  where  the  axisymmetric  solution  did  not  model  well  the  ;• 

flat  bottom  of  the  orbiter.  ■ 
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NOMENCLATURE 

1 

( 

coefficients  of  a partial  differential  equation 

Cf 

skin-friction  coefficient 

s 

specific  heat  at  constant  pressure 

FVSL 

denotes  fully  viscous  shock  layer 

H 

heat-transfer  coefficient, 

H 

stagnation  enthalpy 

H ^ 

ref 

reference  heat-transfer  coefficient 

! 

j 

j = 0 for  two-dimensional  flow 
j = 1 for  axisymmetric  flow 

k 

thermal  conductivity 

1 

L 

length  of  body 

■' 

M 

Mach  number 

M 

molecular  weight  of  mixture  of  gases 

P 

pressure 

Pr 

Prandtl  number,  Cp  y/k 

q 

heat  transfer 

R 

universal  gas  constant 

Re 

Reynolds  number 

r 

body  radius 

St 

Stanton  number,  q/p  U (H„  - H ) 

' 

s 

coordinate  measured  along  body  surface 

T 

temperature 

TVSL 

denotes  thin  viscous  shock  layer 

U,u 

tangential  component  of  velocity 

V 

normal  component  of  velocity 

xi 

w 


X 

y 

a 

6 

Y 

e 

< 

P 

P 


dependent  variable 
coordinate  measured  along  body  axis 
coordinate  measured  normal  to  body 
angle  of  attack 

tt/2  - (Jl 


i 


ratio  of  specific  heats 

1/2 

Reynolds  number  parameter,  (u^gf/p„U^R^) 
surface  curvature 
coefficient  of  viscosity 
density 

angle  between  body  tangent  and  axis 


Subscripts 

sh 

w 

00 
0 

1 
2 


value  behind  the  shock 
wall  value 
frees tream  value 
stagnation  point  value 
conditions  before  the  shock 
conditions  after  the  shock 
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I.  INTRODUCTION 

One  of  the  most  critical  parts  of  any  space  mission  involving  a return  to 
earth  is  the  reentry  into  the  earth's  atmosphere.  Development  of  reusable 
spacecraft,  such  as  the  shuttle,  requires  even  more  knowledge  of  this  important 
phase  of  space  flight.  Reentry  conditions  include  laminar  hypjrsonic  flows 
which  may  be  simulated  by  various  types  of  models.  In  the  approach  most  com- 
monly used,  the  flowfield  is  divided  into  two  parts,  an  inviscid  outer  region 
and  a viscous  boundary-layer  region.  Many  methods  are  available  for  solving 
both  flow  regions.  Examples  of  inviscid  techniques  are  the  methods  of  Inouye, 
Rakich  and  Lomax  (Ref.  1),  Rizzi  and  Inouye  (Ref.  2)  and  Kutler,  Reinhardt  and 
Warming  (Ref.  3).  For  viscous  boundary-layer  solutions,  methods  such  as  that 
of  Blottner  (Ref.  4)  may  be  used.  This  approach  is  most  appropriate  for  super- 
sonic, high  Reynolds  number  conditions  and  encounters  some  problems  when  applied 
to  low  Reynolds  number  flows.  These  problems  include  displacement-thickness  in- 
teraction and  specification  of  the  boundary-layer  edge  conditions.  Problems 
such  as  the  displacement-thickness  interaction  may  be  approximately  treated  by 
using  second-order  boundary-layer  theory  such  as  the  work  of  Lewis  (Ref.  5). 
Specification  of  the  boundary-layer  edge  conditions  is  more  difficult.  For 
supersonic,  high  Reynolds  number  flows  the  boundary  layer  is  thin  compared  to 
the  shock-layer  thickness,  and  the  boundary-layer  edge  conditions  are  considered 
to  be  the  same  as  the  surface  conditions  of  the  inviscid  solution.  For  hyper- 
sonic, low  Reynolds  number  flows,  the  boundary  layer  is  not  thin  and  the  edge 
conditions  must  be  determined  by  another  method  such  as  streamline  tracking  or 
entropy-layer  swallowing. 
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Another  approach  to  solving  hypersonic  blunt  body  flows  has  been  the  use 
of  the  full  Navier-Stokes  equations.  The  elliptic  nature  of  these  equations 
increases  the  complexity  of  the  solution  procedures  and  restricts  their  appli- 
cation in  the  downstream  direction.  Most  of  the  current  methods  have  been  used 
to  solve  the  stagnation  region  and  have  only  been  applied  downstream  a distance 
of  about  one  nose  radius. 

A third  approach  to  the  reentry  problem  is  viscous  shock-layer  flow  in 
which  the  entire  flowfield  from  the  body  to  the  shock  is  treated  in  a uniform 
manner.  Use  of  a viscous  shock-layer  approach  avoids  problems  such  as  stream- 
line tracking  and  displacement-thickness  interaction.  The  equations  are  para- 
bolic in  the  streamwise  direction  so  there  are  no  restrictions  on  obtaining 
solutions  in  the  downstream  direction,  and  at  high  Reynolds  numbers  the  equations  i 

tend  to  the  first-order  boundary-layer  equations. 

The  purpose  of  this  work  has  been  to  validate  a viscous  shock-layer  method  ] 

developed  by  Miner  and  Lewis  (Ref.  6)  by  comparing  solutions  using  that  method  | 

with  experimental  data.  Miner  and  Lewis  followed  the  formulation  of  Davis  i 

(Refs.  7 and  8)  and  Moss  (Ref.  9)  to  develop  a system  of  governing  equations 
for  viscous  shock-layer  flows  of  reacting  gases.  A similar  code  was  developed 
which  used  a perfect  gas  model  and  eliminated  the  complicated  chemical  reactions 
and  species  conservation  equation.  The  perfect  gas  model  was  used  in  the  present 
work  because  it  required  much  less  computing  time  than  did  the  reacting-gas 
code  and  modeled  well  the  wind-tunnel  conditions.  A large  set  of  experimental 
data  was  made  available  by  the  Arnold  Engineering  Development  Center  (AEDC)  where 
tests  had  been  conducted  on  models  of  the  shuttle  orbiter  (Refs.  10-13).  These 
tests  were  conducted  in  two  of  the  tunnels  at  the  von  Karman  Gas  Dynamics  Facil- 

1 

ity  (vKF)  which  could  simulate  conditions  similar  to  those  encountered  during 

reentry.  j 
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Tunnel  B is  a continuous,  closed-circuit,  variable  density  tunnel  which 
can  be  operated  at  stagnation  pressures  from  20  to  300  psia  for  a nominal  Mach 
number  of  6 and  from  50  to  900  psia  for  Mach  8.  The  Reynolds  number  range  was 
0.52  X 10^  to  3.8  X 10^  per  foot.  The  model  was  injected  into  the  tunnel  for 
tests  and  then  retracted  for  cooling  or  changes.  Data  measured  in  Tunnel  B 
at  M = 7.9  included  wall  pressure  distributions,  pitot  pressure  and  tempera- 
ture profiles  (Ref.  10)  and  heat-transfer  distributions  (Ref.  11)  at  angles 
of  attack  from  0 to  45  degrees.  The  second  tunnel  used  was  the  vKF  Hyperveloc- 
ity Wind  Tunnel  (F)  which  is  an  arc-driven  hotshot  wind  tunnel.  This  tunnel  is 
capable  of  Mach  numbers  from  7.5  to  20  and  Reynolds  numbers  per  foot  from  0.75 
10^  to  1.0  X 10^.  The  tests  were  conducted  using  the  nozzle  for  M = 10.1  to 
16.2  and  nitrogen  as  the  test  gas.  Useful  test  times  for  Tunnel  F were  approxi 
mately  100  msec.  Data  taken  in  this  tunnel  were  heat-transfer  distributions 
along  the  body  of  the  model  at  various  angles  of  attack  (Ref.  12).  Calcula- 
tions were  made  using  the  perfect  gas  method,  and  the  computed  results  were 
compared  with  experimental  data  from  these  two  tunnels. 
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In  the  present  work,  the  governing  equations  for  the  viscous  shock-layer 
flows  follow  the  formulation  of  Davis  (Refs.  7 and  8).  The  equations  were  first 
nondimensionalized  by  variables  of  order  one  at  the  body  surface  (corresponding 
to  high  Reynolds  number,  boundary- layer  flows).  The  equations  were  also  non- 
dimensionalized by  variables  of  order  one  in  the  outer  inviscid  flow  (corres- 
ponding to  the  shock  region).  A single  set  of  equations  was  then  obtained  by 
retaining  terms  from  the  equations  in  each  set  to  second  order.  The  resulting 
set  of  shock-layer  equations  was  uniformly  second-order  accurate  in  the  inverse 
Reynolds  number  parameter  e from  the  body  to  the  shock.  Both  longitudinal 
and  transverse  curvature  were  included.  As  given  by  Davis,  the  governing  vis- 
cous shock-layer  equations  were  specialized  for  a perfect  gas  (Ref.  7)  or  a 
binary,  reacting  mixture  of  oxygen  atoms  and  molecules  (Ref.  8).  Miner  and 
Lewis  (Ref.  6)  gave  the  shock-layer  equations  for  a multicomponent  mixture  of 
reacting  gas. 


Governing  Equations 

The  equations  for  shock-layer  flows  of  perfect  gases  are  given  below. 
Continuity  Equation: 


9s 


(r  + y cos  (J))'’  pu 
s-Momentum  Equation: 
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(1  + Ky)  (r  + y cos  (J))"^  pv 


= 0 


1 ^ . 9u  K ^ 1 9P  _ 2 9 

1+Ky  9s  9y  1+Ky  l+<y  9s  " ^ 9y 


1 1 

9u 

ku  j 

H 

9y  ■ 

l+<yj_ 

. 2 
+ E U 


2k  , j cos  4>  1 1 

9u 

KU 

, 1+Ky  r + y cos  0 / ' 

. 9y  ■ 

1+Ky 

(1) 


(2) 


! 

i 


-4- 


I 

I 


! 

I 

■ 


I 

5 

i 

5 

I 

I 

I 

't 


k 

■■t 

I 


y-Momentum  Equation: 


9P  < 2 

3y  1+Ky 


pu  ^ - pv  (FVSL) 
1+icy  ' 3s  3y  ' 


which  becomes 


3P  _ K ,2 
3y  l+<y  P 


(TVSL) 


if  the  thin  shock-layer  approximation  is  made. 


(3a) 


(3b) 


Energy  Equation: 
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Equation  of  State: 
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These  equations  are  in  nondimensional  variables  which  are  defined  by  the 
following  relations; 
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= u /U 
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= V /U 
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For  the  finite-difference  solution  procedure,  it  was  helpful  to  transform 
the  shock-layer  equations.  The  variables  were  normalized  by  the  local  shock 


values  as  follows: 


n = y/y* 


C = s 


u = u/u. 


V = v/v. 


P = P/P 


p = p/p. 


T = T/T, 


P = p/p. 


k = k/k. 


When  written  in  the  transformed  n coordinates,  the  s-momentum  and  energy 
equations  could  be  expressed  in  the  standard  form  for  a parabolic  partial 
differential  equation. 

^ ^ + /\  lli+Au  + A + /\  ^=0  (6) 

3n 

where  W represents  the  dependent  variable.  The  coefficients  A^.  were  functions 
of  the  independent  and  dependent  variables  and  are  given  in  Ref.  6 for  reacting 
gas  mixtures  and  in  Ref.  7 for  perfect  gases. 

Boundary  Conditions 


At  the  body  surface,  the  no-slip  boundary  conditions  were  imposed,  thus 


for  n = 0 


u = 0 


V = 0 


T = T. 


where  T was  either  a specified  constant  or  distribution, 
w 
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At  the  shock,  the  velocity  components  tangent  and  normal  to  the  shock  are 
not  the  same  as  the  components  tangent  and  normal  to  the  body.  The  transforma- 
tions relating  the  shock  components  with  the  body  components  were 


Ush  = i^sh  + 6)  + cos(a  + b) 


Vsh  = -*Jsh  + 6)  + sin(a  + B) 

where  and  are  the  shock  oriented  components  and  B = ir/2  - (j>.  The 
boundary  conditions  at  the  shock  were  determined  by  the  modified  Rankine- 
Hugoniot  relations. 


"sh  = 

a COS  a / ^sin 

. 2 

“ " ^h 

Oi/Sn)^^] 

(10) 

(u^^  - cos  a)^  + 4y  Sin^a/(v+l)^  + |^2/y-1  - 4(y-1 )/(y+1 

]/ 

M^- 

00 

4/(y+1 M^sin^aj 

^ sin  a/2 

sin 

(11) 

Psh 

= ^2  sin^a  - (y- 

•D/yMf]  /( 

Y+1  ) 

(12) 

Psh 

= TP3,/(y-1)  T^, 

(13) 

'^sh 

= -sin  a/p^^ 

(14) 

The  shock  boundary  conditions  for  the  normalized  dependent  variables  at 
n = 1 were  lT  = 7=  ^ = 'p=T=1. 
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Perfect  Gas  Property  Model 
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The  thermodynamic  and  transport  properties  for  the  perfect  gas  were 
calculated  using  standard  relations.  The  temperature  was  computed  by  the 
energy  conservation  equation.  The  density  was  then  calculated  from  the  equation 
of  state  and  the  viscosity  from 

y = (1  + C‘)  (T)^-^(T  + C)  (15) 

where 

C = 198.6  °R/(y-  1)  T„ 

The  transport  properties  were  specified  by  setting  the  Prandtl  number  and 
Lewis  number  to  0.7  and  1.0,  respectively. 

Surface  Conditions 

The  surface  skin-friction  and  heat-transfer  rates  were  given  by  the  skin- 
friction  coefficient  and  Stanton  number.  The  skin-friction  coefficient  in  the 
nondimensional  variables  was  defined  as: 


The  Stanton  number  in  the  nondimensional  variables  was 

St  = q„/(Hg  - (17) 

where 

q„  ■ [k  >T/»y]  „ 

Method  of  Solution 

The  s-momentum  and  energy  equations  were  expressed  in  the  standard  form 
for  a parabolic  partial  differential  equation. 
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+ At  + A, 
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These  equations  were  solved  using  the  algorithm  described  by  Davis  (Ref,  7) 
and  Miner  and  Lewis  (Ref.  6).  The  algorithm  assumes 


w"  = E u"*'  + F„ 
n n 


(19) 


is  valid  through  the  shock  layer  where 


E = 


n + A E , 
n n n- 1 

and 

p °n  - "n  ^-1 

" ■ ^ * "n  En-1 

where  A„,  C„  and  D„  are  as  defined  in  Ref.  6. 
n n n n 


The  boundary  conditions  at  the  wall  were  specified  by  defining  F.j 


and  E-|  such  as 


F,  = W and  E,  = 0 if  W,  = W., 
1 w 1 1 w 


or 

F,  -OandE,  - 1 If  (f  ) ^ =0 

At  the  shock,  the  value  of  W is  set  equal  to  the  shock  value.  Thus  the  solution 
of  the  differential  equation  is  obtained  by  evaluating  E^  and  F^  from  n = 0 
to  n = as  in  Ref.  6,  and  then  using  Eq.  19  to  evaluate  the  dependent 
variable  with  n decreasing  from  to  0. 

Solutions  for  the  continuity  and  n-momentum  equations  were  obtained  by 
integration  with  the  trapezoidal  rule.  The  normal  momentum  equation  is 
rewritten  in  the  transformed  coordinates  so  that  3P/3n  may  be  evaluated  directly 
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as 


9n 


k J'sh 

^sh  ^sh 

p U^- 

"sh 

3V 

Psh  (' 

* k J'sh  ■') 

*^sh 

^ 9n 

^sh  ‘’sh 

^h '' 

i*!-  D IT 

/ 1 7 ^'^sh 

DV  \ 

Psh  ('  ‘ 

'^sh 

r-  U 

n) 

Vsh  DC 

^sh 

Dn  / 

(20) 

The  pressure  derivative  with  respect  to  n is  calculated  and  then  integrated 
to  provide  the  solution  of  the  normal  momentum  equation.  The  continuity 
equation  was  integrated  as  in  Ref.  6 to  yield  the  normal  velocity  profile  and 
the  shock-layer  thickness. 


Solution  Procedure 

At  each  s location  the  shock-layer  equations  were  solved  in  the  following 
order;  energy,  s-momentum,  continuity  and  n-momentum.  The  solution  was 
iterated  until  the  temperature  and  tangential  velocity  converged  at  each 
point  in  the  finite-difference  grid.  After  a converged  solution  was  reached, 
the  profiles  were  used  as  initial  profiles  for  the  next  s location. 

Since  the  equations  depend  upon  dy^^/dc  and  av/8C,  they  possess  an 
elliptic  nature  and  require  global  iterations  to  obtain  an  accurate  solution. 
The  elliptic  effect  in  the  n-momentum  equation  is  resolved  by  considering  TVSL 
flows  for  the  first  global  iteration  and  using  the  calculated  v profiles  for 
subsequent  FVSL  iterations.  The  elliptic  effect  due  to  dy^^/d^  is  resolved 
by  making  a suitable  approximation  for  dy^^/d?  for  the  first  iteration.  If 
the  geometry  was  such  that  Newtonian  flow  was  appropriate,  the  first  iteration 
assumed  dy^^/dc  = 0 and  used  the  calculated  distribution  in  the  following 
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iterations.  If  the  hypersonic  flow  was  non-Newtonian,  an  initial  shock  shape 
was  determined  from  a blunt  body,  method  of  characteristics  procedure  such  as 
that  of  Inouye,  Rakich  and  Lomax  (Ref.  1),  and  the  calculated  dy^^/d^  distri- 
bution was  used  in  subsequent  global  iterations. 


■ 


III.  RESULTS  AND  DISCUSSION 


Solutions  from  the  perfect  gas,  viscous  shock-layer  program  described 
briefly  in  the  preceding  section  were  compared  with  experimental  data  from  the 
AEDC  tests.  Tables  I-IV  list  the  conditions  for  each  of  the  experimental  tests 
used  for  comparison  with  this  numerical  method.  The  data  consist  of  heat- 
transfer  distributions,  wall  pressure  distributions  and  flowfield  surveys. 

Adams  et  al.  (Ref.  13)  describe  the  shuttle  orbiter  windward  streamline  at 

30- deg  angle  of  attack  as  a 31 -deg  hyperboloid.  Figure  1 compares  the  actual 
shuttle  orbiter  geometry  at  'x  = 30-deg  with  the  31-deg  hyperboloid  geometry. 

The  analytical  approximation  models  the  geometry  very  well  to  a distance  of 
x/L  = 0.8  where  the  orbiter  begins  to  curve  upward.  Since  this  31-deg  hyper- 
boloid was  the  only  geometry  available,  it  was  used  as  the  input  geometry  for 
the  shock-layer  program.  Beyond  x/L  = 0.8,  the  geometries  were  not  similar; 
therefore,  comparisons  were  not  made  with  experimental  data  in  this  region. 

To  calculate  the  geometry  of  the  shuttle  at  different  angles  of  attack,  the 

31- deg  hyperboloid  was  rotated  about  the  center  of  a unit  sphere  located  at  the 
nose.  A fifth-deg  polynomial  curve-fit  was  used  to  insure  the  proper  bluntness 
and  curvature  of  the  nose  of  the  rotated  body. 

The  first  data  used  for  comparison  was  the  Tunnel  B body  pressure  distri- 
butions as  given  by  Carter  and  Martindale  (Ref.  10).  Figure  2 shows  the  pres- 
sure distributions  at  various  angles  of  attack  and  a unit  Reynolds  number  of 

4 

2.3  x 10  /cm.  The  agreement  between  the  numerical  prediction  and  the  experi- 
mental data  was  excellent.  At  the  lower  angles  of  attack,  especially  a = 15-deg, 
there  was  an  overexpansion  around  the  nose  of  the  shuttle  and  subsequent  re- 
compression on  the  body.  This  could  not  be  modeled  with  the  assumption  of  a 
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Newtonian  pressure  distribution,  thus  the  numerical  procedure  predicted  an 
expansion  around  the  nose  and  a constant  pressure  along  the  body.  Figure  3 
shows  the  pressure  distribution  for  « = 30  deg  but  at  a higher  unit  Reynolds 

4 

number  of  3.6  x 10  /cm.  Again  the  agreement  with  experimental  data  was  ex- 
cellent. 

Comparisons  with  the  total  temperature  profile  data  of  Martindale  (Ref. 

10)  are  shown  in  Figs.  4-8.  These  profiles  were  taken  at  angles  of  attack 

4 

ranging  from  15-  to  35-deg  and  a nominal  unit  Reynolds  number  of  2 x 10  /cm. 

In  the  experiment,  the  temperature  probe  began  its  survey  at  some  distance 
away  from  the  body  and  traversed  the  shock  layer  by  moving  inward.  During  the 
run,  the  orbiter  model  was  heated  considerably,  especially  near  the  nose,  so 
when  the  probe  neared  the  body,  the  constant  wall  temperature  assumption  was 
no  longer  valid.  Since  only  three  wall  temperature  measurements  were  taken 
along  the  body,  it  was  difficult  to  determine  the  proper  wall  temperature  to 
be  used  in  these  temperature  profile  comparisons.  With  this  consideration  the 
calculated  temperature  profiles  are  in  reasonable  agreement  with  the  experimen- 
tal data.  Figures  4-8  consist  of  two  parts  and  show  the  development  of  the 
profiles  from  x/L  = 0.1  to  x/L  = 0.8.  Figure  4 presents  the  profiles  at 
a = 15-deg  and  shows  a remarkable  difference  between  the  experimental  profiles 
at  x/L  = 0.7  and  at  x/L  = 0.8.  The  numerical  predictions  did  not  show  this 

large  difference.  This  may  be  due  to  an  experimental  problem  or  the  effects  j 

of  the  orbiter  curvature  which  begins  around  x/L  = 0.8.  A higher  unit  Reynolds 

4 

number  solution  is  shown  in  Fig.  7 with  a = 30-deg  and  unit  Re  = 3.7  x 10  /cm. 

The  agreement  for  these  cases  was  as  good  as  those  at  the  lower  Reynolds  numbers. 

Calculations  were  made  for  comparison  with  the  pitot  pressure  profiles 
of  Ref.  10.  Figures  9-13  show  the  pressure  profiles  for  the  same  cases  that 
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were  used  with  the  temperature  profile  comparisons.  Figure  9 presents  the 
profiles  at  a = 15-deg  for  x/L  from  0.1  to  0.8.  The  prediction  of  the  shock 
location  was  considerably  low,  and  the  higher  pressures  were  predicted  in 
the  range  of  x/L  = 0.3  to  x/L  = 0.7.  This  may  be  caused  by  the  Newtonian 
pressure  distribution  used  in  the  first  iteration.  At  low  angles  of  attack 
the  flow  is  non-Newtonian  and  an  overexpansion  and  recompression  occur  as 
seen  in  the  body  pressure  distribution. 

Another  cause  of  the  poor  agreement  at  the  low  angles  of  attack  is 
at  these  conditions  the  axisymmetric  hyperboloid  does  not  model  the  shuttle 
orbiter  geometry  very  well.  Figure  14  compares  two  cross-sections  of  the 
orbiter  with  the  present  model  at  various  angles  of  attack.  Figure  14  shows 
the  orbiter  at  30-deg  and  0-deg  angles  of  attack.  At  30-deg  the  axisymmetric 
hyperboloid  is  a good  model  of  the  orbiter's  flat  bottom.  At  0-deg  the 
radius  is  much  smaller  and  does  not  model  the  flat  bottom  well.  Figures  15 
and  16  show  similar  schematics  at  a = 15-,  25-  and  45-deg.  The  profiles  at 
a = 25-deg  shown  in  Fig.  10  are  in  much  better  agreement  with  the  experimental 
data  than  were  the  profiles  at  a = 15-deg.  At  the  higher  angles  of  attack, 
the  numerical  procedure  begins  to  overpredict  the  shock  location  downstream 
of  the  blunt  nose.  This  may  be  due  to  crossflow  terms  which  have  not  been 
considered  in  this  analysis.  Crossflow  would  cause  a thinner  shock  layer  and 
thus  a lower  experimental  shock  location  than  that  predicted  by  the  present 
viscous  shock-layer  model. 

A number  of  calculations  were  done  using  a two-dimensional  model  to  see 
if  the  two-dimensional  solution  would  give  a better  prediction  of  the  shock 
location  at  low  angles  of  attack.  Flowfield  profiles  were  taken  from  the 
shock-layer  solutions  and  compared  with  inviscid  solutions  using  the  method 
of  Black  and  Lewis  (Ref.  14)  and  with  available  experimental  data.  The  profiles 
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calculated  were  pitot  pressure,  total  temperature,  static  temperature,  tangential 
velocity,  density,  static  pressure  and  Mach  number.  Figures  17-30  show  the 
results  of  both  axisymmetric  and  two-dimensional  solutions  at  a = 15-deg.  The 
agreement  between  the  viscous  and  inviscid  solutions  was  better  with  the  2-D 
model  than  with  the  axisymmetric  model.  The  shock-layer  thickness  predicted  by 
the  2-D  solutions  was  much  greater  than  the  experimental  data.  Even  though  the 
axisymmetric  solution  did  not  model  the  bottom  of  the  orbiter  at  low  angles  of 
attack,  it  still  gave  much  better  predictions  than  the  two-dimensional  model. 
Figures  31-37  compare  the  axisymmetric  viscous  and  inviscid  solutions  for 
a = 30-deg.  The  two-dimensional  model  would  not  produce  a converged  solution 
at  the  higher  angles  of  attack.  The  viscous  solution  is  again  in  reasonable 
agreement  with  the  inviscid  method  of  Black  and  Lewis. 

Figures  38-45  show  heat-transfer  results  from  the  present  method  and  the 

experimental  data  of  the  Tunnel  B convective  heat-transfer  tests.  Figure  38 

shows  the  nondimensional  heat  transfer  at  a ' 0-deg.  The  predicted  values  are 

much  higher  than  the  experimental  data  since  the  present  method  predicted  a 

thinner  shock  layer  at  the  lower  angles  of  attack.  Figures  39-44  compare  the 

heat-transfer  distributions  at  various  angles  of  attack  and  Reynolds  numbers/cm 

4 4 

ranging  from  1.64  x 10  to  9.84  to  10  . The  agreement  with  the  experimental 
data  was  good.  At  the  higher  angles  of  attack,  i.e.  35-  and  45-deg,  the  pre- 
dicted values  were  slightly  lower  since  the  present  method  could  not  model  the 
crossflow.  As  the  Reynolds  number  increased,  the  transition  moved  upstream 
toward  the  nose.  In  Figs.  43  and  44,  the  large  increase  in  the  heat-transfer 
distribution  is  caused  by  the  transition  to  turbulent  flow.  Since  the  present 
model  was  for  laminar  flows  only,  it  could  not  accurately  predict  the  heat- 
transfer  coefficient  in  the  transitional  or  turbulent  region.  Figure  45  shows 
the  upstream  movement  of  the  transition  regime  with  increasing  Reynolds  number 
at  a = 30-deg. 
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Comparisons  with  the  data  of  Tunnel  F are  shown  in  Figs.  46-51.  The 
agreement  with  Tunnel  F data  was  good.  The  present  method  underpredicted  the 
heat-transfer  distribution  along  the  body  because  Tunnel  F used  a nitrogen 
test  gas  which  could  not  be  accurately  modeled  with  the  perfect  gas,  viscous 
shock-layer  program.  Figures  46-51  show  the  predictions  of  the  present  method 

4 

at  various  angles  of  attack  and  Reynolds  number/cm  ranging  from  1.97  x 10  to 

5 

1.57  X 10  . Again  the  movement  of  the  transition  regime  upstream  with  in- 
creasing Reynolds  number  is  apparent  in  these  figures. 

Some  of  the  computing  times  required  for  these  tests  are  given  in  Table  V. 
These  times  were  obtained  on  the  IBM  370/158  system  at  the  Virginia  Polytechnic 
Institute  and  State  University.  The  times  are  for  the  perfect  gas  viscous 
shock-layer  program  used  in  the  present  work  with  a convergence  test  of  one 
percent  on  the  velocity  and  temperature  profiles,  an  initial  step  size  of 
As  = 0.2  and  a maximum  step  size  of  = 1.0. 
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IV.  CONCLUSIONS 


The  results  of  the  present  work  have  shown  that  the  perfect  gas  viscous 
shock-layer  method  may  be  used  to  obtain  accurate  predictions  of  the  blunt- 
body  reentry  vehicle  in  relatively  short  computing  times.  Comparisons  with  data 
from  the  AEDC/vKF  Wind  Tunnel  tests  show  that  the  present  method  makes  good  pre- 
dictions of  the  wall  pressure  distribution,  stagnation  temperature  profile, 
pitot  pressure  profile  and  the  heat-transfer  distribution.  Comparisons  at  the 
low  angles  of  attack  show  that  the  axisymmetric  hyperboloid  does  not  model  well 
the  flattened  bottom  of  the  orbiter  which  causes  a smaller  shock-layer  thickness 
and  higher  heat-transfer  distribution.  At  low  angles  of  attack  the  axisymmetric 
solution  still  gave  better  predictions  of  the  shock  layer  than  a two-dimensional 
solution.  Finally,  at  the  high  angles  of  attack,  solutions  using  the  present 
method  predicted  thicker  shock-layer  thicknesses  and  lower  heat-transfer  dis- 
tributions possibly  because  of  the  neglected  crossflow  terms.  More  work  is 
needed  to  see  if  the  addition  of  the  crossflow  effects  will  improve  the  solu- 
tion at  high  angles  of  attack. 


i 
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TABLE  II 


AEDC  TUNNEL 

B WALL 

PRESSURE 

TESTS 

GROUP 

L 

M 

00 

a 

Re/ft  X 10® 

(ff^) 

■^B 

NO. 

(slugs/ft^) 

(°R) 

(DEG) 

(°R) 

30 

2.272 

98.2 

7.95 

30 

1.110 

600.0 

31 

1.378 

98.9 

7.92 

30 

0.669 

565.0 

32 

1.380 

99.1 

7.92 

35 

0.669 

545.0 

33 

1.383 

98.8 

7.92 

25 

0.672 

545.0 

34 

1.387 

98.3 

7.92 

15 

0.675 

555.0 

TABLE 

III 

AEDC  TUNNEL  B FLOWFIELD  SURVEYS 

GROUP 

NO. 

Poo 

(slugs/ft^) 

T 

on 

(°R) 

M 

00 

a 

(DEG) 

Re/ ft  X 10® 

{ft'h 

^B 

(°R) 

8 

1.377 

99.4 

7.92 

30 

0.685 

605.0 

18 

1.392 

98.8 

7.92 

35 

0.676 

590.0 

46 

1.386 

98.0 

7.92 

25 

0.675 

585.0 

50 

1.368 

99.0 
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AEDC  TUNNEL  F HEAT  TRANSFER  TESTS 
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TABLE  V.  COMPUTING 

TIMES 

CONDITIONS 

NO.  OF 
STATIONS 

NO.  OF 
GLOBAL 
ITERATIONS 

TOTAL 
NO.  OF 
ITERATIONS 

EXECUTION 

TIME 

MIN.-SEC 

TUNNEL  B 

180 

4 

427 

4:30 

TUNNEL  F 


178 


426 


4:09 


Hyperboloid  Geometry  with  the  Windward  Streamline  of  the  Shuttle  Orbiter. 


fp  -0.5334  cm  (0.21  in) 

Pq2*  0.09  atm  (1.323  psia) 

Unit  Re  -2.3  X 10^/cm  (0.69x  lO^/ft) 
L - 57.35  cm  (22.58  in) 

T„  -743.3  °K  (1338.0  ° R) 


TUNNEL  B BODY  PRESSURE  DISTRIBUTION 


=7.95,  = 0.5334  cm  (0.21  in) 

Unit  Re=  3.6 X 10^/cm  (l.lx  lO^/ft) 


T^/Tq  = 0.4199,  = 0. 148  atm  (2. 177  psia) 

L = 57.35  cm  (22.58  in),  Tq  =743.89  °K  (1339.0  °R) 


Figure  3:  Wall  Pressure  Distribution  along  the  Windward  Streamline  of  the 

Shuttle  Arbiter  at  AEDC  Tunnel  B Conditions,  Re/cm  - 3.6  x 10^. 


TOTAL  TEMPERATURE  PROFILES  x/L:  0. 1-0.4 


m 'X 


Fioure  4:  Total  Temperature  Profiles  alonn  the  Windvard  Streamline  of  the  Shuttle  Orbiter  at 

AEOC  Tunnel  B Conditions,  i = 15-dec. 


TOTAL  TEMPERATURE  PROFILES  x/L:  0. 1-0.4 


Figure  5:  Total  Temperature  Profiles  along  the  '.'indward  Streamline  at  25-deg  Angle  of  Attack. 
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PRESSURE  PROFILES 


Continued 


57.35  cm  (22.58  in) 


PITOT  PRESSURE  PROFILES  x/L  0. 1-0.4 


Pressure  Profiles  along  the  Windward  Streamline  at  30-deg  Angle  of  Attack 


Figure  11;  Continued. 
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PITOT  PRESSURE  PROFILES  x/L;  0. 1-0.4 


Figure  13:  Pitot  Pressure  Profiles  along  the  Windward  Streamline  at  35-deg  Angle  of  Attack 
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Cross  Section 


Figure  14:  Cotnparison  of  the  Cross  Sections  of  an  Axisyninetric  Hyperboloid  and  the  Shuttle  Orbiter 

at  a = 30-deg  and  a = 0-deg. 
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SYMMETRIC  INVISCID  AND  VISCOUS  SHOCK  LAYER  TEMPERATURE  PROFILES 


figure  21:  Axisyimetric  Static  Temperature  Profiles  at 


2-0  INVISCIO  AND  VISCOUS  SHOCK  LAYER  TEMPERATURE  PROFILES  AT  a = 15 


•dimensional  Static  Temperatu 
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AXISYMMETRIC  INVISCID  AND  VISCOUS  SHOCK  LAYER  DENSITY  PROFILES 


Figure  25:  Axisymmetric  Density  Profiles  at 


2-D  INVISCID  AND  VISCOUS  SHOCK  LAYER  DENSITY  PROFILES  AT  a = 15 


Figure  26:  Two-dimensional  Density  Profiles  at 


AXISYMMETRIC  INVtSCIO  AND  VISCOUS  SHOCK  LAYER  PRESSURE  PROFILES  AT  a = 15 


Figure  27:  Axisynmetric  Preisure  Profiles  at  a = 15  deg. 


2-D  INVISCID  AND  VISCOUS  SHOCK  LAYER  PRESSURE  PROFILES  AT  a = 15 


Figure  28;  Two-dimensional  Pressure  Profiles  at 


SYMMETRIC  INVISCIO  AND  VISCOUS  SHOCK  LAYER  PITOT  PRESSURE  PROFILES 


Pressure  Profiles  from  an  Axisymnetrlc 


AXISYMMETRIC  INVISCID  AND  VISCOUS  SHOCK  LAYER  TANGENTIAL  VELOCITY  PROFILES 


Tangential  Velocity  Profiles  from  an  Axi symmetric  Solution 


AX  1 SYMMETRIC  INVISCIO  AND  VISCOUS  SHOCK  LAYER  DENSITY  PROFILES 


Figure  35:  Density  Profiles  from  an  Axisymnetric  Solution  at 


AX  I SYMMETRIC  INVISCIO  AND  VISCOUS  SHOCK  LAYER  MACH  NO.  PROFILES  AT  a 


Figure  37;  Mach  Number  Profiles  from  an  Axisymmetric  Solution  at  a = 30-deg. 


■e  43:  Heat-Transfer  D 

Re/cm  = 8.2  x li 


Figure  49:  Tunnel  F Heat-Transfer  Distribution  along  the  Windward 

Streamline,  Re/cm  = 7.55  x 10'^. 
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Figure  51:  Tunnel  F Heat-Transfer  Distribution  along  the  Windward 

Streamline,  Re/cm  = 15.7  x 10^. 
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